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Abstract 

We compute the viscosity spectral function of the dilute Fermi gas for different values of the s-wave 
scattering length a, including the unitarity limit a — > oo. We perform the calculation in kinetic theory 
by studying the response to a non-trivial background metric. We find the expected structure consisting 
of a diffusive peak in the transverse shear channel and a sound peak in the longitudinal channel. At zero 
momentum the width of the diffusive peak is ujq ~ (2e)/(3r/) where e is the energy density and 7] is the 
shear viscosity. At finite momentum the spectral function approaches the collisionless limit and the width 
is of order ujq ~ k(T '/to) 1 / 2 . 



I. INTRODUCTION 



Q3- it 



The dilute Fermi gas at unitarity is a strongly correlated scale invariant quantum fluid 
provides a beautiful model system for testing ideas from many different areas of physics, including 
condensed matter physics, nuclear and particle physics, and quantum gravity An important 
property of the Fermi gas at unitarity is nearly perfect fluidity J3L Experiments with trapped 
Fermi gases in the vicinity of a Feshbach resonance indicate that [4|t9| 

'/ 



s 



<0.5h/k B , (1.1] 



where r\ is the shear viscosity and s is the entropy density. This value is smaller than the viscosity to 



entropy density of any known fluid with the possible exception of the quark gluon plasma 



10I-I12] 



Equ. (11.11) approaches a bound rj/s = h/i^Anks) that has been proposed based on the AdS/CFT 
correspondence jisl. \l^. 

The Kubo formula relates the shear viscosity to the zero momentum, zero frequency limit of the 
retarded stress tensor correlation function. In this work we will study the frequency and momentum 
dependence of this correlation function. We will focus on the regime of high temperature and 
small frequency and momentum where the methods of kinetic theory are applicable. Our study is 
motivated by several considerations. The first observation is that the frequency and momentum 
dependence of the retarded correlator encodes the time and distance scales over which the viscous 
contribution to the stress tensor relaxes to the Navier-Stokes expression. These scales are important 
or understanding how hydrodynamics breaks down in the dilute corona of a trapped Fermi gas 
3, [3] • Kinetic theory predicts that in the dilute limit the shear viscosity is independent of density. 
This leads to a difficulty in the analysis of experiments involving scaling flows in which the velocity 
field is linear in the coordinates. The Navier-Stokes equation predicts that in this case the viscous 
stress tensor is a function of time, but constant in space. As a consequence the dissipative force, 
which is proportional to gradients of the stresses, vanishes and the dissipated heat, which is governed 
by the volume integral of the square of the stress tensor, is infinite. This problem can be resolved if 
a finite relaxation time for the dissipative stresses is taken into account. In this case the dissipative 
stress tensor vanishes in the dilute corona, the dissipated heat is finite, and there is a viscous force 
that counteracts the expansion of the system. 

A second motivation is to understand the physics behind nearly perfect fluidity in the unitary 
Fermi, in particular the question whether the unitary gas can be described in terms of quasi- 
particles. In kinetic theory quasi-particles manifest themselves in terms of a peak in the spectral 
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function, but the quasi-particle peak may disappear in strong coupling. The spectral function is 
accessible through the behavior of the euclidean correlation function which can be computed in 
imaginary time Monte Carlo simulations. A calculation of the spectral function in kinetic theory 
serves as a useful default model for numerical calculations of the euclidean correlator. 

This work extends and complements several recent studies of the transport properties of the 
dilute Fermi gas at unitarity. Calculations of transport coefficients in the context of kinetic theory 



can be found in 
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-|20l]. and a quasi-particle model is discussed in 2l|. A diagrammatic calculation 



of the shear viscosity using the Kubo formula is described in |22| . This work also studies the spectral 



function. A sum rule for the spectra 
on dilute Bose gases can be found in 
to non-relativistic systems are described in 
in QCD is studied in (2TH31 1 . 



function of the stress tensor was proved in 



Related work 



24( , and first attempts to extend the AdS / CFT correspondence 
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26] . Finally, the spectral function of the stress tensor 



II. PRELIMINARIES 



The retarded correlation function of the stress tensor is defined by 

G% kl (u;,k) = -i fdtf dxe^ lk - x e(t)([n ij (t,x),n H (o,o)]). 



(2.1) 

A careful definition of the stress tensor U a "(t, x) in terms of the Hamiltonian of the dilute Fermi 



gas was given m 



22 



32| . In the following, we only need the definition of the stress tensor in terms 



of hydrodynamic and kinetic variables, see Sec. II III and HVl The spectral function is defined by 



P 



ijkl 



u,k) = -2ImG ijkl (uj,\t) 



(2.2) 



In the following, we will focus on the longitudinal and transverse shear channels p xyxy , p xzxz anc l 
pzzzz where rhe ^-direction is the direction defined by the external momentum k. We will compute 
the retarded correlation function using linear response theory. The external field that couples 
to the stress tensor is the metric gij(t, x). A formalism for coupling an interacting theory of 
non-relativistic particles to an external metric in a way that exhibits a non-relativistic version of 
general coordinate invariance was recently developed in 33[ . The response to a small perturbation 



1.1 



hij(t,x) around the flat metric is given by 
5U ij 



(2.3) 



'1.1 



where II|? is the expectation value of the stress tensor in equilibrium. 
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III. HYDRODYNAMICS 



At very low frequency and momentum the retarded correlation function is governed by 
hydrodynamics. The hydrodynamic limit provides an important constraint for the kinetic 
description of the system. In this section we will derive these constraints in the case of correlations 
functions of the stress tensor. We note that hydrodynamics breaks down for frequencies u > (r]/n)T 
and momenta k > (r)/n) (mT) 1//2 , where n is the density and T is the temperature. In the unitary 
limit and in the vicinity of T c we have rj/n ~ 1 and the range of validity of hydrodynamics is very 
large. In the weak coupling limit 77 ~ (mT) 1 / 2 /a 2 , where a is the scattering length, and the range 
of validity of hydrodynamics is much smaller. 

In order to study linear response theory we consider hydrodynamics in a curved background 



(7ij(x, t). Non-relativistic fluid dynamics in a non-trivial metric was studied in 34|. The continuity 
equation and the equations of momentum and energy conservation are given by 

-±-d t (y/gp) + V, (p<) = 0, 
^(V^O+V.IF = 0, 

^jft (y/gps) + v, (W - £#r) = ^ , (3.1) 

where p = mn is the mass density, v l is the fluid velocity, and g = det(gij). The covariant derivative 
of a vector field is given by ViVj = diVj — F^Vk, where • is the Christoffel symbol associated with 
gij. The stress tensor is 

IF = pv i v j + Pg lJ - a ij , (3.2) 
where the viscous term, aij, is given by 

<Tij = V (V^- + VjVi + g i3 ^ + U - IjA 9ij (v k v k + . (3.3) 

The dissipative function R is given by 

2R = 2 ( ViUi + VjVi ~ % 9i i Vkvk + fa ~ i^'f) 

+C( V ^ + ^) +f9 l Td l T, (3.4) 

where k is the thermal conductivity. We are interested in small deviations from equilibrium. For 
this purpose we write g^ = 5ij + hij, and p = p + 5p, P = P + 5P etc. and linearize the equations 



in hij, Vi, 8P, 5s, 5p and ST. After going to Fourier space we find 

^-up h + iu8p + p (-ik)v z = 0, 

iup v i + {-ik)8P8 iz + ika iz = , 
nk 2 

icup 5s + —ST = , (3.5) 

where h = Tr(hij) and we have chosen the direction of k to be the z-axis. The fluctuations in 
thermodynamic variables are related by thermodynamic identities. We will consider 8p and 8s to 
be the independent variable and use 

dT, T 

Op ,s Cy 

where c s is the speed of sound and cy is specific heat at constant volume. 

We can now solve for the hydrodynamic variables in terms of the external field and compute 
Ily from Eq. ( 13.21) . We then determine Gr using Eq. (12.31) . We are particularly interested in the 
transverse shear correlators G^ xy and G^ xz , as well as the longitudinal correlator G Z R ZZZ . We find 
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Gxzxz _ ITjUJ . . 

R — ., 2 f I \ ' \ 6 -°> 

uj -ik 2 (r]/p ) 
G zzzz = uJ 2 po(c 2 s +iua) 
" 2 k 2 - \uk 2 " 



dP 


\s S P + 


dP, 


dp 


~ds~\p 


dT 


l>+ 


dT, 


dp 


ds \p 



UJ 2 — C 2 K' 



We note that G x ^ xz has a diffusive pole, where r]/po is the momentum diffusion coefficient, and 
G Z B ZZZ has a sound pole where a is the sound attenuation coefficient, 

4 KC 2 k 2 ( 1 1 \ . 
aA> = ^ + C+— V • (3-10) 



3 ' LO 2 \cy Cp 

Finally, we can take the imaginary parts and compute the spectral functions. In the Kubo limit 
k = 0, u — y we find 

lim J± = U + C, (3.12) 

w,fc->o leu 3 

where T, L labels the transverse and longitudinal components of the spectral functions, pt = 

pxyxy ^ pxzxz an( j ^ _ ^zzzz _ 



IV. KINETIC THEORY 



In this section we calculate the spectral function using kinetic theory. The basic object in kinetic 
theory is the distribution function /(x, p, t) where f = j\ = fi- The distribution function satisfies 
the Boltzmann equation. The kinetic equation in a curved background can be found by starting 
from the Boltzmann equation in general relativity 
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^ - r > v !0 x ' p) = m ' (4 - i} 



where i,j,k are three-dimensional indices and u, a,/3 are four-dimensional indices. In the non- 
relativistic limit p° ~ m, Tqq ~ 0, and rL ~ \g lk gkj [33(- We get 

(I + iw - + r '^P) w) m x ' " = c|/l • (4 ' 2) 

We consider small deviations from equilibrium and write / = fo + Sf with /o(p) = 
/ (p l p J ^y/(2mT)). We also write = Sij + and linearize in and Sf. We get 

(dt + P -d^ Sf + ^f^PYk = C[Sf] . (4.3) 

We will solve the linearized Boltzmann equation by making an ansatz for 5f. We go to Fourier 
space and restrict ourselves to quadratic terms in p and write 

£/(£*;, k,p) =w/o(l-/ )— = J - — r— r-, 4.4 

2mi w — v p • k + le 



where £tl = (tl(w, k) and 

1 L 1 

Using this ansatz the LHS of the linearized Boltzmann equation becomes 



L ^ = {(& + U/'V/'i + (& + 1)pV^} • (4-6) 

The RHS of the Boltzmann equation involves the linearized collision integral 

RHS = J dr 234 r(l, 2; 3, 4)^ [^( Pl ) + ^(p 2 ) - <^(p 3 ) - ^(p 4 )] (fr^ + ^g) , (4.7) 

where we have defined <7 tJ '(q) = g* J (q; a;, k) = q l q> /{uj — \ q ■ k + ie) and we have labeled the 
momenta such that pi = p. We have also defined the phase space measure dr 2 34 = dridr 2 dr 2 
with dr, = d 3 pi/(27r) 3 and the transition rate r(l,2;3,4) = w(l, 2; 3, 4)£>(1, 2; 3, 4). The factor 
.D(l,2; 3, 4) contains the distribution functions 

£>(1, 2; 3, 4) = / (pi)/o(p 2 ) (1 - /o(p 3 )) (1 - /o(p 4 )) , (4.8) 



and to (1, 2; 3, 4) is the collision probability 

w{l, 2; 3, 4)oT 3 4 = v rel dfi' . (4.9) 

Here, v re i = \\\ — v 2 | is the relative velocity, dQ' is the solid angle along 2q' = p 3 — p 4 and 

(da)/(dQ') is the differential cross section. We also define 2q = pi — p2 where |q| = |q'|. We will 

consider a dilute Fermi gas in which the scattering amplitude is completely characterized by the 

s-wave scattering length a. The differential cross section is 

da 1 a 2 , > 

dfi 7 ~ 47 f + a?q 2 ' ' 

In order to solve for £t,l we will take moments of the linearized Boltzmann equation. The quadratic 
moments of the LHS are given by 

L a b = I ^ iUjf f m - f0) p« P W £ k ¥&,L + 1) , (4.H) 

where L ab = L ab (u, k). The moments of the RHS are given by 
^ = 2~T \ dV ^ r ^ 2 ^ A )P^ b d9^ 

m J T,L 

dT tot r(l, 2; 3, 4) [p*p\ + p%p\ - p%p\ - p a iP \] «^'( Pl ) £ &, L h\ 



2mT 

2u f d 3 pd 3 q a 2 q 

^202 \q~q- - q~q~ ) g ' [Pi) 
with D 2 ^ 2 = D(l, 2; 3, 4) and R ab = R ab (co, k). 



T,L 



^fjwF d!V Tvk* ° 2 " 2 (?V ~ « V) 9 ' J(PI) £ • (412) 



Transverse channel 



We can now determine £t,l by solving L ab = R ab . We begin by considering h xy = h yx ^ with 
all other components of vanishing. In this case h = and only £t contributes. We find 

LXV = / dr ^ ((T + 1)pW ^ 

= ^ + 1)^ f dpp 6 /o(l - /o) - ^±12^ (4.13) 
307r 2 mT 7 FF J0K J0J 15 y K J 



and 



2u£ T h xy f d 3 pd 3 q a 2 q .. , 

2u£ T h xy f a 2 q p 2 q 2 q x q y q x q y _ uj£, T h x , 



Ti 3 mT 



x v f a a a ar a 1 n p 2 q 2 q x q y q x q v _ u£ T h xy 
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Setting L xy = R xy gives 

The superscript (xy) indicates that we have obtained £t by projecting on the U xy channel. In the 
limit k — y our ansatz Eq. (14.41) is complete, and there is no dependence on the channel. At finite 
k, however, additional tensor structures in Eq. (14. 4p are possible and there is some dependence of 
^L t T on the channel. 

In kinetic theory the stress tensor is given by 

where the factor 2 is the spin degeneracy. Using the ansatz Eq. (14. 4p we get 

mT J (2nf Joy J0J mu -p z k + ie ^ T y y V 1 

We can now extract the retarded correlation function and the spectral function using Eq. (I2.3P and 
Eq. (|22). We find 

(fy*v(u,k) = 2lm(^ v \u,k))Re(j xy (u,k)) + 2Re(^ v \u,k))lm(j xy (u,k)) . (4.18) 

In general the functions C xy and J xy have to computed numerically. The limit k — > can be 
studied analytically. We may consider two cases. In the free case the collision integral C xy vanishes 
and £t = — 1- The spectral function is then given by p xyxy (u, k) = — 21m J xy (co, k). The function 
J xy (u, k) can be computed analytically for all values of u and k. The collisionless spectral function 
is 

p xyxy (u,k) 1 (2mT) 2 ( x ( mu 2 \ 



2u 167r k 

where Li2(x) is a polylogarithm and (f = exp(— fx/T) is the fugacity. In the limit k — > we find 

o xyxy (u k) e 

lim / J"'** = y5{u), (4.20) 

k->0 2bJ 6 

where e is the energy density 

d 3 p p 2 



At finite momentum the delta function is spread out over the regime to < \k\(T/m) 1 / 2 . In the 
collisional case we find that in the limit k — y the collision integral is real and scales as C xy (u, k — y 
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0) ~ 1/uj. This means that the spectral function is dominated by the lm(^ y \uj, k))Re( J xy (oo, k)) 
term. We also find that in this limit Ke(J xy ) ~ I xy /(15m). This implies that 

V( U ) = lim pXyX ^ k) = -&> MC xy (u,0)) (4 22) 

The zero frequency limit is 

r? = lim = rr = — yL- I , 4.23 

^ 15m.Re(C xy (u,0)) 32^F | 1/(3mTa2) a ^ Q 

where in the last step we have evaluated the integrals in the high temperature limit. This result 
agrees with the hydrodynamic limit given in Eq. (13. 71} and the known formula for the shear viscosity 



of a dilute Fermi gas in kinetic theory 



161 ] . We can now write the frequency dependence as 



1 + CU Z T Z 



where r = I xy /(uTle(C X y)) = 15mrj/I xy is a relaxation time. In the high temperature limit r agrees 



with the result obtained in 



15j. We can also evaluate the sum rule 



dur](uj) = ^ = ^ = -. (4.25) 
ir J 2r 30m, 3 



We note that the sum rule agrees with the integral of the free spectra 



function given in Eq. (I4.20p . 
23| if we take account the fact 



The sum rule also agrees with the general sum rule obtained in 
that kinetic theory does not reproduce the short distance contribution 77(0;) ~ C/^mu, where C 
is Tan's contact [37 1. 

The spectral function in the xz channel can be computed in the same fashion. The spectral 
function p xzxz (u, k) is given by Eq. (I4.18P with the obvious replacement (xy) — > (xz). In the Kubo 
limit we find 

lim pXZX ^ k) = V (4.26) 

cj,k->0 2u 



as expected from hydrodynamics, see Eq. (13. 8j) . Numerical results for the spectral function in the xy 
and xz channel are shown in Fig. [TJ We consider two different scattering length, one corresponding 
to weak coupling, kpa = 0.1, and one close to unitarity, kpa = 100. We plot the spectral function 
as a function of oj/T for several values of k = kj\j2mT and T/Tp. 

We first discuss p xyxy (u, k). At zero momentum the spectral function normalized to 2770; is a 
Lorentzian with a width proportional to I/77. This implies that the peak in the spectral function 
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FIG. 1: Stress tensor spectral functions in the xy and channels. We plot the spectral function p(u),k) 
normalized to the value in the Kubo limit, p(uj — > 0, 0) = 2t]uj, as a function of uj/T for different values of 
k = k/ (2mT) 1 / 2 and T/Tp. The panels on the left show the result for kpa = 0.1, which corresponds to 
weak coupling, and the panels on the right show the spectral function close to unitarity, for kpa = 100. 

is narrow in weak coupling, but very broad in strong coupling. As the momentum k grows the 
width of the peak increases and it becomes less sensitive on the strength of the coupling. In this 
regime the spectral function is close to the free spectral function. We also find numerically that 
the spectral function satisfies the sum rule Eq. (14.251) even at non-zero momentum k. In the xz 
channel the normalized spectral function p xzxz (u, k)/(2rya;) vanishes in the limit u — > for any 
non-zero value of k. For k = the normalized spectral function approaches unity as oj — > 0. This is 
in agreement with the hydrodynamic prediction Eq. f)3.8p . We find that the spectral function has 
a peak which, like the peak in the xy channel, corresponds to a diffusive mode. The width of the 
mode is controlled by k except in the limit of small momenta where the width is controlled by the 
shear viscosity. 
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B. Longitudinal case 



In this section we study the longitudinal spectral function p zzzz (u, k). The calculation is more 
involved because both £y and £x contribute. We consider a perturbation h zz ^ with all other 
hy = 0. In this case hfj = and h T zz = - = -^f- diag(l, 1, -2). The LHS of the 

Boltzmann equations is L zz = Lf? + L Z L Z with 



2mT 



(4.27) 



The collision integrals can be written as 



2u£ T f d 3 pd 3 q a 2 q , , •■ T 

^2f J -7CT d ^^ 2 l , _ a _J gY-|/y*)g ,3 (Pi)hS 



1 + erg 



/ , J _ J _ J _. JO „ « 2 <? ^ 2 „ 4 (/3 2 -|)(p 2 (l-3a 2 )-8 M a/3 + 4g 2 (l-3/3 2 )) 



127r 3 mT 
uj£ T h 



dp dq da d/3 D 2 ^2 



1 + a 2 g 2 



p g 



2moj — pka — 2qk(3 + ie 



T 5 



m 2 T 7 (2vr 



1 + erg 



37r 3 mT 



iii ion A 2 4 (/3 2 -|)(p 2 /4 + pga/3 + g 2 ) 
dp dg da d/J D 2 ^2 — — ^ P q 



1 + a 2 g 2 2mu — p£;a — 2qk(3 + ie 



(4.28) 



(4.29) 



We can now solve for £t,l- We find 



TZZ 

1 T.L 



The zz component of the stress tensor is given by 



5U Z 



2lo 



d 3 p £M ^PVPY ^hjj + ^hfj 



(27T) 



/o(l - /o) 



m 2mT lo — v D • k + ie 



dp da 



/o(l - /o)p 6 « 2 



4 22) (3« 2 -l)+e 



(zz) 



12ir 2 mT J r moo — pka + ie L" 
= cuh zz (^ ) ( W ,k)J^( W ,k)+^ ) (a;,k)Jf(a;,k)) . 



(4.30) 



(4.31) 



We can now extract the retarded correlator and the spectral function using Eq. (I2.3P and Eq. (12.21) . 

We get 

p zzzz {u,k) 



2cu 



Im ^\u, k)J£>, k) + & z \u, k) Jf (u, k) 



(4.32) 



The result again simplifies in the limit k — > 0. In the collisional case we find that the function J is 



10 




FIG. 2: Viscosity spectral function in the zz channel. We plot the spectral function p(cj, k) normalized to 
the value in the Kubo limit, p(u — > 0, 0) = ^rjco, as a function of u/T for different values of k = fc/(2mT) 1 / 2 
and T/Tp. 

dominated by its real part. Also, one can show that the k — > limit of Re(C£ z ) = and therefore 
Im(£i zz) ) = 0. We find 

lim P ZZZZ ^ k K li m I m (^ ) (a;,k))Re(j^(c,k)) = - f^^^fe^ (4-33) 
k^o 2u k^o [l z T z ) 2 + Re(Cf (w,0)) 2 v ' 

The transverse and longitudinal integrals that we have done are obviously related to each other. In 
the limit k — > we have 

2 4 4.2 

I? = —Ixy J T z (u, 0) = - J xy (uj, 0) = — I xy C T Z = —C xy . (4.34) 

This shows that as k — > the ansatz Eq. (I4.4p is complete and £y (oj) = ^^(oo). We can also 
show that 

Um = i (435) 

o;,k^0 2W 3 V ^ 

Comparing this result with the hydrodynamic prediction (I3.9D we find that the bulk viscosity C, 



vanishes. This follows from scale invariance at unitarity [34J. The fact that we also find Q = 
away from unitarity is related to the fact that we only take into account elastic 2 <H> 2 collisions. 
Numerical results for p zzzz (u,k) normalized to the hydrodynamic limit |w?7 are shown in Fig. [2j 
We observe that the spectral function shows a sound peak at u = c s k. At the largest temperature 
considered, T/Tp = 3, the speed of sound differs from the asymptotic high temperature value 
(? s = (5T)/(3m) by about 10%. The width of the sound peak is close to that of the free spectral 
function. This is at variance with hydrodynamics which predicts that the sound peak is narrow in 
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strong coupling, and in the limit of small momentum. At zero momentum hydrodynamics predicts 
that p{uj)/uj is constant. Higher order effects will turn this part of the spectrum into a broad peak 
with a width inversely proportional to the relaxation time. This implies that at small momentum 
we expect a narrow peak with a width proportional to (k/c s )(r)/p) superimposed on a broad peak 
with a width proportional e/(Trj). Our ansatz correctly reproduces the k — > limit as well as the 
existence of a sound peak, but it appears to be too simple to reproduce the two-peak structure at 
small k. 



V. SUMMARY AND DISCUSSION 

We have computed the viscosity spectral function using kinetic theory in a non-trivial background 
metric. We find that a simple quadratic ansatz for the non-equilibrium distribution function 5f 
satisfies the main constraints on the spectral function, in particular the existence of a diffusive peak 
and a sound peak (although it misses the interplay between these two features in the longitudinal 
channel), the correct Kubo limit, and the viscosity sum rule. Kinetic theory misses the non- 
analytic high frequency tail of the spectral function 77(0;) ~ C/^mu. We find that the width of 
the diffusive peak at zero momentum is inversely proportional to the viscosity. This implies that 
the relaxation time for the shear stress is r ~ 3r)/ (2e), in agreement with [l5|. We also find that at 
large momentum the spectral function is close to the free spectral function, and the width of the 
transport peak is of order kiT/m) 1 ^ 2 . 

There are at least two important questions that we have not addressed in this paper. The first is a 
systematic matching to second order hydrodynamics. For this purpose it is important to understand 
the structure of second order hydrodynamics in the conformal limit along the lines of the analysis 



in the relativistic case performed in [38[. The second is a detailed study of the imaginary time 
correlation functions, with an emphasis on finding the optimal strategy for extracting transport 
coefficients from Quantum Monte Carlo calculations. The main difficulty is that the imaginary 
time correlation function is mainly sensitive to the spectral weight of the transport peak, and that 
the sum rule implies that the weight is independent of the shear viscosity. This problem can be 



addressed by utilizing the information contained in the momentum dependence of t 



re correlation 



function, and by considering not only the shear channel but also the sound channel 27J, l30 |. 
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